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Abstract 

Plasmodium of Physarum polycephalum is a single cell visible by unaided eye. During 
its foraging behaviour the cell spans spatially distributed sources of nutrients with 
a protoplasmic network. Geometrical structure of the protoplasmic networks allows 
the plasmodium to optimize transfer of nutrients between remote parts of its body, 
to distributively sense its environment, and make a decentralized decision about 
further routes of migration. We consider the ten most populated urban areas in 
United Kingdom and study what would be an optimal layout of transport links 
between these urban areas from the "Plasmodium's point of view". We represent 
geographical locations of urban areas by oat flakes, inoculate the plasmodium in 
Greater London area and analyse the Plasmodium's foraging behaviour. We simulate 
the behaviour of the plasmodium using a particle collective which responds to the 
environmental conditions to construct and minimise transport networks. Results of 
our scoping experiments show that during its colonization of the experimental space 
the plasmodium forms a protoplasmic network isomorphic to a network of major 
motorways except the motorway linking England with Scotland. We also imitate 
the reaction of transport network to disastrous events and show how the transport 
network can be reconfigured during natural or artificial cataclysms. The results of 
the present research lay a basis for future science of bio-inspired urban and road 
planning. 
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1 Introduction 



Plasmodium of Physarum polycephalum\^\ is a single cell with many diploid 
nuclei. The plasmodium feeds on microbial creatures and microscopic food 
particles. When colonising its habitat the plasmodium tries to optimize the 
network of its protoplasmic protrusions and tubes to span the maximal number 
of nutrient sources, and to minimize costs of transportation and intra-cellular 
communication complexity. At the same time the plasmodium positions its 
body away from the domains of chemo-repellents, low humidity and high il- 
lumination. Due to its highly sophisticated adaptive behaviour, distributed 
sensing and actuating, and decentralized decision making, the plasmodium is 
considered to be amongst the most prospective experimental prototypes of 
biological amorphous computers [T8|19|I20|I21|I22] . 

Decision-making routines inside the Plasmodium's body are executed by in- 
teracting bio-chemical and excitation waves [18], redistribution of electrical 
charges on Plasmodium's membrane pQ and spatio-temporal dynamics of me- 
chanical waves [TS]. The Plasmodium's 'intelligence' is based on the interac- 
tions between propagating patterns. Therefore the plasmodium can be consid- 
ered as an adaptive intelligent react ion- diffusion [3] or excitable [T1I7] medium 
encapsulated in elastic growing membrane — a novel type of biological com- 
puting architecture at the interface of reaction-diffusion chemical processors [2] 
and amorphous computers [9]. Experimental proofs of P. polycephalum compu- 
tational abilities include approximation of shortest path [21J and hierarchies 
of planar proximity graphs [6], computation of plane tessellations [28J, im- 
plementation of primitive memory [26], execution of basic logical computing 
schemes [31J, control of robot navigation [32], and natural implementation of 
spatial logic and process algebra [27] . 

Nature-inspired computing paradigms and experimental implementations have 
already been successfully applied to the approximation of a minimal- distance 
path between two given points of a two-dimensional space or a road network. 
The shortest-path problem is solved in experimental reaction-diffusion chem- 
ical systems [2J, gas-discharge analog systems |25j, spatially extended crystal- 
lization systems [S] , and using computer and mathematical models of collective 
insects [TU] and P. polycephalum [29] . Intriguing analogies between mycelian 
fungi networks and road networks are discussed in [13] . Previously [5] we have 
evaluated a road-modeling potential of P. polycephalum, however no conclusive 
results were presented there. In the present paper we made a full-scale — and 
as you will see further successful — 'attack' on the following Physarum-ro&d- 
building problem. Given major cities represented by oat flakes and plasmod- 
ium of P. polycephalum inoculated in one of the cities, will the plasmodium 
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develop a protoplasmic network connecting oat flake that matches the net- 
work of motorways connecting the cities? We show that for UK motorways, 
connecting ten most populous urban areas, the answer is largely 'y es ' with 
some interesting exclusions. 

The paper is structured as follows. We describe experimental and simulation 
setup in Sect. [2} In Sect. [3] we discuss experimental and simulation results 
on developing transport links between major urban areas by plasmodium. 
Scoping experiments on plasmodium imitating reconfiguration of road links 
and mass migration caused by disastrous accidents in selected urban areas 
are reviewed in Sect. |4j Section [5] suggest lines of further experiments and 
highlights possible applications of plasmodium-computing. 



2 Methods 

Plasmodium of P. polycephalum is cultivated in plastic container, on paper 
kitchen towels sprinkled with distilled water and fed with oat flake^" 2 "] For 
experiments we use 120 x 120 mm polyestyrene square and 90mm diameter 
round Petri dishes. We use either 2% agar gel (Select agar, Sigma Aldrich) or 
a moistured filter paper as a non-nutrient growth substrate. Agar plates and 
filter papers are cut in a shape of United Kingdom island. 



We consider the ten most populous urban areas in United Kingdom (Fig. PT 



— Greater London, Bristol, Sheffield, Nottingham, Liverpool, Tyneside, Greater 
Glasgow, West Yorkshire, Greater Manchester, and West Midlands — as per 
2001 CensuQ see details and boundaries of the areas in 1241. The areas are 
projected onto gel or filter paper and oat flakes (of size approximately match- 
ing size and shape of the areas) are placed in the positions of the urban areas 
(Fig. (TJd). At the be ginning of each experiment plasmodium is inoculated in 
the centre of Greater London Urban Area. 

The Petri dishes with substrate and plasmodium are kept in darkness, at tem- 
perature 22-25 C°, except for observation and image recording. Periodically 
the dishes are scanned using Epson Perfection 4490 scanner. Large-scale dis- 
asters in urban areas were imitated by salt crystals^ being placed in centres 
of the urban areas. 

Scanned images of dishes are enhanced for higher visibility, saturation in- 
creased to 204, and contrast to 40. To ease readability of experimental images 

2 Asda's Smart Price Porridge Oats 

3 Office for National Statistics, General Register Office for Scotland and Northern 
Ireland Statistics and Research Agency 

4 Saxa corse grain sea salt, RHM Foods, CW10 0HD, UK 
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Fig. 1. Experimental basics, (a) Schematic map of ten most populous urban areas, 
shown in black, (b) Typical experimental setup: urban areas are represented by oat 
flakes, plasmodium is inoculated in London, the plasmodium spans oat flakes by 
protoplasmic transport network. Map of urban areas is adapted from |24j . 

we provide complementary binary version of each image. The binarization is 
done as follows. Each pixel of a color image is assigned black color if red R 
and green G components of its RGB color exceed some specified thresholds, 
R > 9 R , G > 9 G and blue component B does not exceed some threshold 
value B < 6b] otherwise, the pixel is assigned white color (exact values of the 
thresholds are indicated in the figure captions as 9 = 9c, 9b)- 

The computational approximation of Physarum is provided by the particle 
collective approach introduced in [H] where a population of very simple mo- 
bile particles with chemotaxis-like sensory behaviour were used to construct 
and minimise spatially represented emergent transport networks in a diffusive 
environment. The particle approximation corresponds to a particle approxi- 
mation of LALI (Local Activation Long-range Inhibition) reaction-diffusion 
pattern formation processes [H] and exhibits a complex range of pattern- 
ing by varying particle sensory parameters [15]. We assume that each parti- 
cle in the collective represents a hypothetical unit of Physarum plasmodium 
gel/sol interaction which includes the effect of chemoattractant gradients on 
the plasmodium membrane (sensory behaviour) and the flow of protoplas- 
mic sol within the plasmodium (motor behaviour). The summation of particle 
positions corresponds to a static snapshot of network structure whilst the col- 
lective movement of the particles in the network corresponds to protoplasmic 
flow within the network. 

Although the model is very simple in its assumptions and implementation it is 



4 



capable of reproducing some of the spontaneous network formation, network 
foraging, oscillatory behaviour, bi-directional shuttle streaming, and network 
adaptation seen in Physarum using only simple, local functionality to gen- 
erate the emergent behaviours. Details of the particle morphology, sensory 
and motor behavioural algorithms can be found in [T7] and in this paper we 
use an extension of the basic model (without utilising oscillatory behaviour) 
to include plasmodium growth and adaptation (growth and shrinkage of the 
collective). 

Growth and adaptation of the particle model population is currently imple- 
mented using a simple method based upon local measures of space availabil- 
ity (growth) and overcrowding (adaptation, or shrinkage, by population re- 
duction). This is undoubtedly a gross simplification of the complex factors 
involved in growth and adaptation of the real organism (such as metabolic in- 
fluences, nutrient concentration, waste concentration, slime capsule coverage, 
bacterial contamination etc.). However the simplification renders the popula- 
tion growth and adaptation more computationally tractable and the specific 
parameters governing growth and shrinkage are at least loosely based upon 
real environmental constraints. Growth and shrinkage states are iterated sep- 
arately for each particle and the results for each particle are indicated by 
tagging Boolean values to the particles. The growth and shrinkage tests were 
executed every two scheduler steps and the method employed is specified as 
follows. 

If there are 1 to 10 particles in a 9 x 9 neighbouhood of a particle, and the 
particle has moved forwards successfully, the particle attempts to divide into 
two if there is an empty location in the immediate neighbourhood surrounding 
the particle. If there are to 24 particles in a 5 x 5 neigbourhood of a particle 
the particle survives, otherwise it is annihilated. 

To replicate the experimental spatial configuration the environment is repre- 
sented by a colour coded image. Urban areas are indicated by black areas, 
uninhabitable coastal boundary of the UK mainland is represented by the 
light grey border region and the background (empty substrate or nutrient-rich 
substrate, depending on the experiment) is shown as white. We simulate three 
environmental conditions. The first is a nutrient rich-substrate with back- 
ground nutrients in addition to oat flakes, corresponding to Physarum growth 
on oatmeal agar. The second condition is a high-concentration nutrient source 
on a non-nutritious substrate. This corresponds to growth on plain agar. The 
third condition is foraging behaviours which approximates dendritic foraging 
growth on damp filter paper. 

The three different environmental conditions are approximated by adjusting 
simulated chemoattractant projection strengths and diffusion distances in the 
particle model. Oat flake positions were represented by projection of chemoat- 
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tractant to the diffusion map at the locations of urban areas. The projection 
weight (altering the strength of chemoattractant projection, higher weight val- 
ues result in more chemoattractant projection to the environment) applied to 
the urban areas was 0.01 in the nutrient-rich and foraging experiments, and 
was set to 0.5 for the high-concentration experiments. Diffusion damping (af- 
fecting the distance of chemoattractant gradient diffusion, lower values result 
in longer diffusion distances) was set to 0.1 for nutrient-rich substrate and 
foraging experiments, and to 0.01 for the high concentration experiments. 

The diffusion kernel size was 7x7 pixels for all experiments. We assume that 
diffusion of chemoattractant from any urban area is suppressed when the area 
is covered by particles. The suppression is implemented by checking each pixel 
within an urban area and reducing the projection value, or concentration of 
chemoattractants, by multiplying it by 0.001 if there is a particle within a 
3x3 neighbourhood surrounding the pixel. 

Contamination of an urban area with salt was implemented by deleting any 
particles which entered the contaminated area. Inoculation of a small initial 
particle population (5 particles) was initiated at random locations within the 
London urban area. 

Particle sensor offset was 5 pixels and angle of rotation set to 45 degrees. 
Particle sensor angle was 30 degrees in the foraging experiments (the relatively 
narrow sensor angle stimulates branching transport networks) and 45 degrees 
in all other experiments. 

The emergent transport networks are represented in the results by the config- 
uration of the particle collective, shown as a spatial map of particle positions. 
Where appropriate, superpositions of previous particle occupancy history are 
shown in grey, and darker shading indicates more frequent transport. Video 
recordings of the dynamical growth and adaptation of the particle model can 
be found at http://uncomp.uwe.ac.uk/jeff/urban/urban.htm, 



3 Development of transport links 



Typically, being placed in the centre of London urban area, plasmodium con- 
sumes some nutrients from its nearest (London) oat flake and starts prop- 
agating North, North- West or West (Fig. |2^ig). Birmingham (Fig. |2^g) and 
Bristol (Fig. (3^d) are the usual candidates which are spanned by London- 
originated plasmodium. When urban area areas concentrated in Midlands are 
colonised by plasmodium and linked by protoplasmic tubes the plasmodium 
heads North towards the Tyneside urban area (Fig. |2j)h). After taking on 
Tyneside the plasmodium propagates North, cross Scottish boundaries and 
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(a) t =12 h (b) t =23 h (c) t =34 h 




(j) t =47 h (k) t =69 h {£) t =80 h 

Fig. 2. Typical Plasmodium development: (a)-(f) scanned image of experimental 
Petri dish; (g)-(l) binarized images, 6 = (200, 200, 150) 
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(a) t =63 h (b) t =63 h (c) t =47 h 




(d) t =47 h (e) t =53 h (f) t =53 h 




(g) t =46 h (h) t =46 h 



Fig. 3. Example of protoplasmic networks which occurred in experiments. Results 
of four independent experiments are shown in (ab), (cd), (ef) and (gh), respectively. 
Time passed after inoculation of plasmodium in London is shown in sub-figure 
caption, (a), (c), (e), (g) are scanned image of experimental Petri dishes; (b), (d), 
(f), (h) are binarized images. Thresholds of binarization are = (170, 170, 100) for 
(b) and (d), and 9 = (130, 130, 100) for (f) and (h). 

finally reaches Glasgow urban area (Fig. (2^i). Then plasmodium continues 
colonization of the substrate until all urban areas (sources of nutrients) are 
colonized (Fig. gl-f and Fig. 

Example of plasmodium networks connecting urban areas are shown in Fig. |3j 
The figures demonstrate that, in general, the structure of the network does 
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(a) t =8 h 



(b) t =31 h 




Fig. 4. Examples of reconfigurations of protoplasmic network: (a) and (b) scanned 
images, (c) and (d) binarised images, = (140, 160, 100) 

not depend significantly on size and shape of the substrate but mainly on 
the configuration of sources of nutrients: 90 mm round Petri dish (Fig. [3^b), 
120 mm rectangular dishes fully covered with agar gel (Fig. [3fc— f ) and shape 
of UK island cut of agar gel plate in 120 mm rectangular Petri dish (Fig. (3j;h). 
Noticeably, the plasmodium does not stop its foraging activity even when all 
sources of nutrients are colonised. It propagates away from the 'designated' 
area (Fig. [3^-f ) unless stopped by an unfriendly substrate (like the bottom of 
a plastic Petri dish not covered by gel, Fig. [3^h). 

The plasmodium does not always keep all sources of nutrients spanned by 
its protoplasmic tubes. Sometimes some tubes cab be abandoned during col- 
onization. Thus in Fig. [4] we see that at the beginning of its development 
plasmodium links London and Nottingham (Fig. |4^lc) . When urban areas in 
Midlands are colonised and linked by protoplasmic tubes to Tyneside and 
Glasgow areas, the plasmodium abandons its tube connecting Greater Lon- 
don and Nottigham urban areas (Fig. [^d). 

The particle approximation results varied depending on the nutrient and dif- 
fusion conditions of the simulated environment. Fig. [5] (a) illustrates radial 
expansive growth in a nutrient-rich substrate environment, approximating the 
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growth on oatmeal agar, for example. The background nutrients (very pale 
grey on areas not reached by population) are assumed to be consumed imme- 
diately on contact whereas the the simulated oat flakes (urban areas) remain 
a persistent source of nutrients. The expansion of the population continues 
until all the substrate area has been convered and then the particle popula- 
tion automatically adapts and shrinks the transport network to encompass 
and connect the urban areas. Fig. [5] (b) is in an environment with a nutrient 
poor substrate but high concentration of chemoattractants emitted from the 
urban areas in highly diffusive conditions. The population expands to locate 
the nearest source of food and then continues until all food sources are located 
and the network between them is minimised. Although this behaviour appears 
to be computationally efficient, it does not reflect the actual reconfiguration 
behaviour seen in Physarum and this is more accurately reflected in Fig. [5] 
(c) where a narrow sensor angle of 30 is used to imitate foraging, branching 
behaviour. The discovery of the food-source nodes (not shown) is not as pre- 
dictable or direct when nutrient conditions are poor and the transport network 
undergoes constant reconfiguration whilst keeping the food sources covered. 

Even a limited number of examples (e.g. Figs. [3]and|4j) demonstrate that Plas- 
modium of P. polycephalum is a very dynamic system, in which the morphol- 
ogy is continuously changing and in which spatio-temporal dynamics rarely 
reaches a fixed stable point (unless humidity decreases and Plasmodium forms 
a sclerotium). There is no such thing as a stationary configuration of proto- 
plasmic network, therefore when extracting a generalised graph of transport 
links from plasmodium experiments we can only build a Physarum-graph P 
as follows. 

Let U be a set of ten most populous urban areas, S={Si, ■ ■ ■ , Sk} is a set 
of series Si, i — 1, . . . , k (k is a number of experiments), of scanned images 
of plasmodium networks, Si = (sj, . . . , s™ 1 ^. For any two areas a and b from 
U the weight of edge (ab) is calculated as follows: w(ab) = J2si&sX(Si,a,b), 
where x{Si, a,b) = 1 if there is at least one snapshot s e 5; which shows 
a protoplasmic tube connecting a and b. We do not take into account exact 
configuration of the protoplasmic tubes but merely their existence. Each pro- 
toplasmic tube is counted just once for any particular experiment. We will also 
consider sub-graphs P a of P, a = 5, 10, 12 defined as follows: for d,i 6 U: 
(ab) G P a if w(ab) > a. 

Physarum-graphs extracted from 25 laboratory experiments are shown in 
Fig. [6j maximum edge weight is 22. The graph becomes planar when we re- 
move edges with weights below 6 (Fig. |6p). The graphs is acyclic, or a tree, 
when only edges appearing in over 40% of experiments are shown (Fig. |6]i). If 
we increase cut-off value to 14 the graph becomes disconnected, and the node 
corresponding to the Bristol urban area becomes isolated. 
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(c) t = 2340, 7736, 11780, 13068, 22552 and 29412 sched- 
uler steps 

Fig. 5. Particle approximations of emergent transport networks with different nu- 
trient concentrations, diffusion properties and foraging behaviour, (a) Evolution of 
radial expansive growth front and subsequent network adaptation on a nutrient-rich 
substrate, (b) evolution of network with a propagating growth front in a highly dif- 
fusive environment, (c) snapshots showing dynamical reconfiguring network when 
foraging behaviour is used. 



11 



(a) 



(b) 




Fig. 6. Physarum-gr&phs for various values of edge weights: (a) all edges of Physarum 
graph are shown, thickness of each edge is proportional to the edge's weight, (a)-(d) 
only edges with weights exceeding 5 (b), 10 (c) and 12 (d) are shown. 

Let us check if there is any correspondence between transport links built by 
Physarum and man-build motorways. We construct the motorway graph M 
as follow^ Let U be a set of ten most populous urban areas, for any two 
areas a and b from U, the nodes a and b are connected by an edge (ab) if 
there is a motorway starting in vicinity of a and passing in vicinity of b and 
not passing in vicinity of any other urban G U. The motorway graph 



5 Graph M is extracted from motorway network as shown in maps . google . com 
and www . openstreetmap . org 
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(a) 

Fig. 7. Graph M of man-made motorway network connecting ten most populous 
urban areas. 

M is shown in Fig. 7 By comparing M (Fig. [7]) and Physarum graphs P and 
its subgraphs (Fig. 6) we found that. 



Finding 1 Motorway graph M is a sub-graph of Phys&rum- graph P. 

This shows that — in princple — the distributed 'logic' underpinning Plas- 
modium's decision-making routines corresponds to human logic behind road- 
planning decisions. However graph P is not planar thus poses no practical 
importance. Motorway link M6/M74 connecting Great Manchester and Great 
Glasgow urban area is represented by plasmodium only in 3 of 25 experi- 
ments. The corresponding edge does not appear in graphs P5, P10 and P12 
(Fig. [6]3 d) . Motorway M4, linking Greater London and Bristol urban areas, 
is represented by Physarum only in 20% of experiments, graph P5 in Fig. [6}d. 

Finding 2 Physarum polycephalum satisfactory approximates motorway net- 
work linking the ten most populous urban areas in United Kingdom except the 
motorway link M6/M74 connecting Greater Manchester and Greater Glasgow 
urban areas. 



Is there a rationale behind Plasmodium's behaviour? Let us have a look at two 
most popular planar proximity graphs, Relative Neighbourhood graph [3U] 
RNG (Fig. §1) and Gabriel graph [TT1I25] GG (Fig. constructed over 
nodes corresponding to centres of urban areas. Points a and b are connected 
by an edge in RNG if no other point c is closer to a and b than dist(a, b) [30] . 
Points a and b are connected by edge in GG is disc with diameter dist(a, b) 
centered in middle of the segment ab is empty [TTf2"3"] . The graphs are related 
as RNG C GG [30.23.12j. Both graphs are imporant in spatial analysis and 
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(a) (b) 

Fig. 8. Two proximity constructed on urban areas U: (a) relative neighbourhood 
graph RNG, (b) Gabriel graph GG. 

statistics, particularly GG which was invented specially to analyse and simu- 
late geographical distribution of biological populations [TT]. Our experiments 
shown that: 

Finding 3 Pi C GG and P 12 = RNG. 

Moreover P12 is a minimum spanning tree MST over U. We know that 
MST C RNG [30J but in the particular case of urban areas' configurations 
U we even have MST(U) = RNG(U). 

With regards to the relationship between motorway graph and proximity 
graphs we see that M is neither sub- nor super-graph of RNG and GG. 
To transform M to RNG one needs to remove two edges from and relocate 
one edge in M, while to transform M to GG it no edges should be removed 
from but three edges added to M. The edge connecting Tyneside and Greater 
Glasgow is present in RNG and GG but absent in M. 

Finding 4 Experiments with plasmodium of Physarum polycephalum show 
that motorway M6/M74 is not optimally positioned and should be rerouted 
from Newcastle to Glasgow. Alternatively M6/M74 may remain intact but new 
motorway Newcastle- Glasgow must be built. 

The evolution and minimisation of transport networks by the particle approx- 
imation in diffusive, nutrient-rich environments closely approximates certain 
variants of the Physarum-graphs, most notably P^- Another similarity be- 
tween Physarum-gr&phs and the emergent transport networks is the common 
absence of the Bristol-London link and the Manchester- Glasgow link. The con- 
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(a) (b) (c) 

Fig. 9. Superpositions of the history of emergent transport networks produced by 
the particle populations. Darker shading indicates more frequently visited — more 
stable — paths, (a) Superposition of network evolution on a nutrient-rich substrate, 
(b) Superposition of evolution in a highly diffusive environment, (c) Superposition 
of evolution of a dynamical reconfiguring network when foraging behaviour is used. 
Note that (a) and (c) show regions (for example in lower left corner, Cornwall) 
where part of the exploring growth front was trapped within the narrow confines of 
the map boundary. 



nectivity of the transport networks is also affected by the nutrient and diffusion 
strength — both the nutrient-rich substrate example and high-concentration 
example show an approximation of a tree structure similar to the RNG. With 
the highly foraging example (Fig. [5^) the connectivity is more varied and dy- 
namic but there is still the tendency to avoid the Manchester-Glasgow and 
Bristol-London links. 



By recording the paths of the particles during the evolution of the emergent 
transport networks is is possible to spatially approximate the variation seen 
in the Physarum-gr&phs. Fig. [9] shows the effect of superimposing the tra- 
jectories of all particles during evolution as a single image. The most visited 
paths are shown as darker colours and the overall trajectory as a cloud-like 
mass around the nodes. Fig. [9|a), formed in by the adaptation of a radially 
expanded population shows remnants of paths formed during the expansion 
process. Fig. |9](b), formed by migratory extension shows a more space efficient 
search (corresponding somewhat to Physarum-graph P\o). The highly foraging 
superposition Fig. |9]^c) shows similarities to Physarum-gr&ph. P 5 and although 
there is significant rearrangement of the network (indicated by the fuzzy na- 
ture of the superposition), the central 'core' of the network (connecting the 
areas of West Midlands, Nottingham, Sheffield, West Yorkshire, Manchester 
and Liverpool) maintains its connectivity. This suggests that foraging is more 
significant in areas away from closely clustered sources of food, even when the 
cluster is far from the initial inoculation site. 
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(a) t =32 h 



(b) t =44 h 



(c) t =59 h 




(d) t =32 h (e) t =44 h (f) t =59 h 



Fig. 10. Plasmodium's response to disaster in, and subsequent contamination of 
West Yorkshire area, filter paper substrate: (a)-(c) scanned images, (d)-(f ) binarized 
images, = (200,200,20). Grain of salt placed (marked by circle in (a)) in West 
Yorkshire urban area at 32 h of Plasmodium development. Images (b)-(c) and (e)-(f) 
show the Plasmodium's response to increased concentration of sodium chloride. 

4 Imitating disasters 



When making experiments it is difficult to resist an impulse to imitate a large- 
scale disaster leading to contamination of one of the urban areas spreading 
to the surrounding areas. A disaster is imitated by placing a grain of salt 
in a substrate's loci corresponding to the urban area. Diffusion of sodium 
chloride in the substrate imitates progressive contamination of surrounding 
areas making them temporarily uninhabitable. 



Here we consider three examples: disaster in and contamination of West York- 
shire urban area (Fig. 



10, wet filter paper substrate, and Fig. 



substrate) and Tyneside urban area (Fig. 12, agar gel substrate) 
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agar gel 



In experiment shown in Fig. 10 a plasmodium network spanning most urban 
areas is formed 32 h after inoculation of Greater London area with Plasmod- 



ium. We place a salt crystal in Leeds (West Yorkshire urban area) (Fig. 10 id). 
In response to increased concentration of sodium chloride the plasmodium 
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(a) t =37 h (b) t =51 h (c) t =85 h 




(d) t =37 h (e) t =51 h (f) t =85 h 



Fig. 11. Plasmodium's response to disaster in, and subsequent contamination of 
West Yorkshire area, agar gel substrate: (a)-(c) scanned images, (d)-(f) binarized 
images, = (200,200,20); a grain of salt placed (marked by circle in (a)) in West 
Yorkshire urban area at 37 h of Plasmodium development. Images (b)-(c) and (e)-(f) 
shows the Plasmodium's reaction toward increased concentration of sodium chloride. 



abandons West Yorkshire area and disconnects the area from neighbouring Ty- 
neside, Greater Manchester and Sheffield urban areas (Fig. flObe) . At the same 
time the plasmodium starts mass-exploration of Scotland and North Wales, 
and restores previously abandoned transport link with London (Fig. [lOpe). 
At around 27-30 h after the 'disaster' in West Yorkshire the plasmodium com- 
pletes its evacuation from Midlands and North England and regroups itself in 
North Scotland and South England (Fig. [lofcf). 



Due to - possibly - a lower rate of sodium chloride diffusion in agar gel (com- 
paring to wet filter paper), the Plasmodium's response to imitated contam- 
ination of agar gel is less dramatic. As in the previous experiment, we wait 
until plasmodium forms a well-established protoplasmic network and then 



place a salt crystal in Leeds (Fig. 11 id). The plasmodium temporarily breaks 
all transport links leading to the contaminated zone (West Yorkshire), and 
increases exploratory activity in Wales (even developing pronounced proto- 
plasmic routes in West Wales) and South- West England (Fig. |llpe). One and 
a half day after contamination of West Yorkshire the plasmodium restores 
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(a) t =32 h 



(b) t =41 h 



(c) t =52 h 






(e) t =32 h 



(f) t =41 h 




52 h 



(h) t =64 h 



Fig. 12. Plasmodium's response to disaster in, and subsequent contamination of 
Tyneside area, agar gel substrate, (a)-(d) scanned images, (e)-(h) binarized images, 
= (200, 200, 20); a grain of salt placed (marked by circle in (a)) in Tyneside urban 
area at 32 h of plasmodium development. Images (b)-(c) and (e)-(f) shows the 
Plasmodium's reaction toward increased concentration of sodium chloride. 

transport links with Leeds (Fig. |TTfcf) and decrease its activity in Wales and 
South- West. 



In the example shown in Fig. 12 we strike urban area with contaminant before 
any transport link leading to the area is established. Plasmodium inoculated 
in Greater London area spans urban areas in Midlands with a protoplasmic 



network (Fig. 12 le). When the Plasmodium's active zone approaches Tyneside 
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area we place salt crystal in Newcastle (Fig. |12pe). In response to contami- 
nation the plasmodium abandons its attempt to colonize North England and 
Scotland and guides its foraging activity towards the West and South, and 9 h 



after contamination of Tyneside the plasmodium reaches Bristol (Fig. 12 Df). 



The plasmodium increases its exploration of Wales and South- West England 



(Fig. 12 zg) and, in 32 h after disaster strikes Tyneside, restores transport 



network between Midlands and London, this time via Bristol (Fig. 12 zg 



Finding 5 Contamination of a single urban area stimulates exploration of 
uncolonized areas and leads to restoration of previously abandoned transport 
links. 

To approximate contamination using the particle model we introduced simu- 
lated salt to the West Yorkshire urban area into a pre-existing foraging trans- 



port network (Fig. 13 (a)). The effect of particle exposure to contaminated 
areas resulted in the removal of those particles from the environment. The con- 
taminated region increased in radius and resulted in disconnection of contam- 
inated nodes from the network. Although this resulted in increased foraging in 
nearby areas, the results were not appreciably different from the general forag- 
ing behaviour of the particle approximation and certainly not as dramatic as 
observed in the organism itself. When the contamination was removed the for- 



aging behaviour resulted in reconnection of the transport network. In Fig. 13 
(b) we simulated the growth of a plasmodium from the London urban area 
in a nutrient-rich substrate where the West Yorkshire urban area was already 
contaminated by salt. The contaminated area 'killed' any particles entering 
the affected region and the growth wave by-passed the region. At t = 1289 the 
contamination was removed and the growth front began to extend into the 
West Yorkshire area, reconnecting the network. After the network was fully 
connected, minimisation of the network occurred in the usual manner. 



5 Discussion 



We undertook scoping laboratory and simulation experiments on bio-inspired 
road planning. We represented the ten most populated UK urban areas by 
source of nutrients, inoculated plasmodium of Physarum polycephalum in one 
of the areas and analysed dynamics of colonizing areas by the plasmodium. 
We performed the same spatial representation of urban areas in a spatially 
represented particle approximation of Physarum under three different envi- 
ronmental conditions: nutrient-rich background substrate, high quality nutri- 
ents in a strongly diffusive environment, and a poor quality environment with 
foraging behaviour. We studied space-time dynamics of spanning the urban 
areas by the Plasmodium's network of protoplasmic tubes and demonstrated 
that the plasmodium transport network sufficiently well matches the topology 
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(a) t = 2449, 2884, 3189, 3655, 3855 and 4372 scheduler steps 




(b) t = 117, 199, 321, 443, 1190, 1539, 2812 and 15000 scheduler steps 



Fig. 13. Effect of disastrous contamination in the West Yorkshire urban area after 
and prior to network formation, (a) Foraging behaviour resulted in connection of 
nodes. Simulation of contamination by Salt with increasing radius. At t = 4228 
contamination is removed and network reconnects with previously contaminated 
nodes, (b) Growth of simulated plasmodium in a nutrient-rich substrate where West 
Yorkshire urban area is already contaminated by salt (radius 37 pixels). At t =379 
the growth front by-passes contaminated area. At t =1289 the contaminated area 
shrinks and the previously contaminated nodes are added to the network. 

of the existing man-made motorway networks with two exclusions: Motorway 
M4 (Bristol-London) rarely occurs in plasmodium networks, and the route 
M6/M74 is entirely absent. Very similar network patterns were observed in 
the particle approximation which also succeeded in minimising the transport 
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network configuration in nutrient-rich environments. As a 'by-product' of the 
experiments we provided an insight into bio-inspired response to disastrous 
contamination of urban areas. Two main components of the response are ex- 
ploring non-colonised territories and restoration of abandoned transport links. 
The approximation of disastrous contamination in the particle model reflected 
the temporary destruction of the transport networks, their reconfiguration af- 
ter the contamination had passed, but did not show evidence of dramatic shifts 
in foraging behaviour when exposed to the contaminant. We believe that the 
lack of dramatic response of the particle collective is because the oscillatory 
behaviour seen in Physarum was not included in this instance of the particle 
collective. We aim to integrate growth behaviours with oscillatory behaviours 
in future works. 

Our experiments did not take terrain into account. This may explain the 
'anomalous' situation with plasmodium not imitating route M6/M74 but de- 
veloping a transport link directly from Newcastle to Glasgow. Terrain based 
experiments on Physarum road planning may form a subject for further stud- 
ies. 
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